Theoretical Study of the Defects and Doping in Tuning the Electrocatalytic Activity of Graphene for CO2 Reduction.

The application of graphene-based catalysts in the electrocatalytic CO2 reduction reaction (ECO2RR) for mitigating the greenhouse effect and energy shortage is a growing trend. The unique and extraordinary properties of graphene-based catalysts, such as low cost, high electrical conductivity, structural tunability, and environmental friendliness, have rendered them promising materials in this area. By doping heteroatoms or artificially inducing defects in graphene, its catalytic performance can be effectively improved. In this work, the mechanisms underlying the CO2 reduction reaction on 10 graphene-based catalysts were systematically studied. N/B/O-codoped graphene with a single-atom vacancy defect showed the best performance and substantial improvement in catalytic activity compared with pristine graphene. The specific roles of the doped elements, including B, N, and O, as well as the defects, are discussed in detail. By analysing the geometric and electronic structures of the catalysts, we showed how the doped heteroatoms and defects influence the catalytic reaction process and synergistically promoted the catalytic efficiency of graphene.

Theoretical Study of the Catalytic Mechanism of the Cu-Only Superoxide Dismutase.

The catalytic mechanisms for the wild-type and the mutated Cu-only superoxide dismutase were studied using the hybrid density functional B3LYP and a quantum chemical cluster approach. Optimal protonation states of the active site were examined for each stage of the catalytic cycle. For both the reductive and the oxidative half-reactions, the arrival of the substrate O2•- was found to be accompanied by a charge-compensating H+ with exergonicities of -15.4 kcal·mol and -4.7 kcal·mol, respectively. The second-sphere Glu-110 and first-sphere His-93 were suggested to be the transient protonation site for the reductive and the oxidative half-reactions, respectively, which collaborates with the hydrogen bonding water chain to position the substrate near the redox-active copper center. For the reductive half-reaction, the rate-limiting step was found to be the inner-sphere electron transfer from the partially coordinated O2•- to CuII with a barrier of 8.1 kcal·mol. The formed O2 is released from the active site with an exergonicity of -14.9 kcal·mol. For the oxidative half-reaction, the inner-sphere electron transfer from CuI to the partially coordinated O2•- was found to be accompanied by the proton transfer from the protonated His-93 and barrierless. The rate-limiting step was found to be the second proton transfer from the protonated Glu-110 to HO2- with a barrier of 7.3 kcal·mol. The barriers are reasonably consistent with experimental activities, and a proton-transfer rate-limiting step in the oxidative half-reaction could explain the experimentally observed pH-dependence. For the E110Q CuSOD, Asp-113 was suggested to be likely to serve as the transient protonation site in the reductive half-reaction. The rate-limiting barriers were found to be 8.0 and 8.6 kcal·mol, respectively, which could explain the slightly lower performance of E110X mutants. The results were found to be stable, with respect to the percentage of exact exchange in B3LYP.

Theoretical investigation on the reaction mechanism of UTP cyclohydrolase.

Nucleoside triphosphate cyclohydrolase (UrcA) is a critical enzyme of the uracil catabolism pathway that catalyses the two-step hydrolysis of uridine triphosphate (UTP). Although the recently resolved X-ray structure of UrcA in complex with substrate analogue dUTP provided insights into the structural characteristics of the enzyme, the detailed catalytic mechanism, including how the reaction intermediate accomplishes conformational conversion in the active centre, remains unclear. In this study, extensive DFT calculations and MD simulations were performed to investigate the catalytic reaction process of UrcA. This study shows that the first hydrolytic reactions in UrcA follow a three-step mechanism, while the second hydrolytic reaction follows a two-step mechanism. Glu392 plays a critical role in deprotonating the lytic water in both hydrolytic reactions. The rate-limiting step of the first hydrolytic reaction lies in the cleavage of the uracil ring, in which an extraneous water molecule bridges the proton transfer from C6-OH to N1 to enable the reaction to go through a six-membered transition state with relatively low steric tension. In the second hydrolytic reaction, Glu392 abstracts protons from the lytic water and directly transfers them to the nitrogen atom of the cleaved C4-N3 bond so that the hydrolytic reaction is no longer rate-limited by the C-N bond cleavage step. MD simulations show that the reaction intermediate experiences spontaneous conformation overturn in the active site of UrcA under the assistance of the hydrogen bond interaction from Tyr307 to place its C4-N3 bond alongside the Zn2+ centre of the enzyme to trigger the second hydrolytic reaction.

Reorienting Mechanism of Harderoheme in Coproheme Decarboxylase-A Computational Study.

Coproheme decarboxylase (ChdC) is an important enzyme in the coproporphyrin-dependent pathway (CPD) of Gram-positive bacteria that decarboxylates coproheme on two propionates at position 2 and position 4 sequentially to generate heme b by using H2O2 as an oxidant. This work focused on the ChdC from Geobacillus stearothermophilus (GsChdC) to elucidate the mechanism of its sequential two-step decarboxylation of coproheme. The models of GsChdC in a complex with substrate and reaction intermediate were built to investigate the reorienting mechanism of harderoheme. Targeted molecular dynamics simulations on these models validated that harderoheme is able to rotate in the active site of GsChdC with a 19.06-kcal·mol-1 energy barrier after the first step of decarboxylation to bring the propionate at position 4 in proximity of Tyr145 to continue the second decarboxylation step. The harderoheme rotation mechanism is confirmed to be much easier than the release-rebinding mechanism. In the active site of GsChdC, Trp157 and Trp198 comprise a "gate" construction to regulate the clockwise rotation of the harderoheme. Lys149 plays a critical role in the rotation mechanism, which not only keeps the Trp157-Trp198 "gate" from being closed but also guides the propionate at position 4 through the gap between Trp157 and Trp198 through a salt bridge interaction.

Insights into the aggregation mechanism of RNA recognition motif domains in TDP-43: a theoretical exploration.

The transactive response DNA-binding protein 43 (TDP-43) is associated with several diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) due to pathogenic aggregations. In this work, we examined the dimer, tetramer and hexamer models built from the RRM domains of TDP-43 using molecular dynamics simulations in combination with the protein-protein docking. Our results showed that the formations of the dimer models are mainly achieved by the interactions of the RRM1 domains. The parallel ß-sheet layers between the RRM1 domains provide most of the binding sites in these oligomer models, and thus play an important role in the aggregation process. The approaching of the parallel ß-sheet layers from small oligomer models gradually expand to large ones through the allosteric communication between the α1/α2 helices of the RRM1 domains, which maintains the binding affinities and interactions in the larger oligomer models. Using the repeatable-superimposing method based on the tetramer models, we proposed a new aggregation mechanism of RRM domains in TDP-43, which could well characterize the formation of the large aggregation models with the repeated, helical and rope-like structures. These new insights help to understand the amyloid-like aggregation phenomena of TDP-43 protein in ALS and FTLD diseases.

Theoretical Studies on the Binding Mode and Reaction Mechanism of TLP Hydrolase kpHIUH.

In this work, we have investigated the binding conformations of the substrate in the active site of 5-HIU hydrolase kpHIUH and its catalytic hydrolysis mechanism. Docking calculations revealed that the substrate adopts a conformation in the active site with its molecular plane laying parallel to the binding interface of the protein dimer of kpHIUH, in which His7 and His92 are located adjacent to the hydrolysis site C6 and have hydrogen bond interactions with the lytic water. Based on this binding conformation, density functional theory calculations indicated that the optimal catalytic mechanism consists of two stages: (1) the lytic water molecule is deprotonated by His92 and carries out nucleophilic attack on C6=O of 5-HIU, resulting in an oxyanion intermediate; (2) by accepting a proton transferred from His92, C6-N5 bond is cleaved to completes the catalytic cycle. The roles of His7, His92, Ser108 and Arg49 in the catalytic reaction were revealed and discussed in detail.

Investigation of the Detailed AMPylated Reaction Mechanism for the Huntingtin Yeast-Interacting Protein E Enzyme HYPE.

AMPylation is a prevalent posttranslational modification that involves the addition of adenosine monophosphate (AMP) to proteins. Exactly how Huntingtin-associated yeast-interacting protein E (HYPE), as the first human protein, is involved in the transformation of the AMP moiety to its substrate target protein (the endoplasmic reticulum chaperone binding to immunoglobulin protein (BiP)) is still an open question. Additionally, a conserved glutamine plays a vital key role in the AMPylation reaction in most filamentation processes induced by the cAMP (Fic) protein. In the present work, the detailed catalytic AMPylation mechanisms in HYPE were determined based on the density functional theory (DFT) method. Molecular dynamics (MD) simulations were further used to investigate the exact role of the inhibitory glutamate. The metal center, Mg2+, in HYPE has been examined in various coordination configurations, including 4-coordrinated, 5-coordinated and 6-coordinated. DFT calculations revealed that the transformation of the AMP moiety of HYPE with BiP followed a sequential pathway. The model with a 4-coordinated metal center had a barrier of 14.7 kcal/mol, which was consistent with the experimental value and lower than the 38.7 kcal/mol barrier of the model with a 6-coordinated metal center and the 31.1 kcal/mol barrier of the model with a 5-coordinated metal center. Furthermore, DFT results indicated that Thr518 residue oxygen directly attacks the phosphorus, while the His363 residue acts as H-bond acceptor. At the same time, an MD study indicated that Glu234 played an inhibitory role in the α-inhibition helix by regulating the hydrogen bond interaction between Arg374 and the Pγ of the ATP molecule. The revealed sequential pathway and the inhibitory role of Glu234 in HYPE were inspirational for understanding the catalytic and inhibitory mechanisms of Fic-mediated AMP transfer, paving the way for further studies on the physiological role of Fic enzymes.

N-Nitrosation Mechanism Catalyzed by Non-heme Iron-Containing Enzyme SznF Involving Intramolecular Oxidative Rearrangement.

The non-heme iron-dependent enzyme SznF catalyzes a critical N-nitrosation step during the N-nitrosourea pharmacophore biosynthesis in streptozotocin. The intramolecular oxidative rearrangement process is known to proceed at the FeII-containing active site in the cupin domain of SznF, but its mechanism has not been elucidated to date. In this study, based on the density functional theory calculations, a unique mechanism was proposed for the N-nitrosation reaction catalyzed by SznF in which a four-electron oxidation process is accomplished through a series of complicated electron transferring between the iron center and substrate to bypass the high-valent FeIV═O species. In the catalytic reaction pathway, the O2 binds to the iron center and attacks on the substrate to form the peroxo bridge intermediate by obtaining two electrons from the substrate exclusively. Then, instead of cleaving the peroxo bridge, the Cε-Nω bond of the substrate is homolytically cleaved first to form a carbocation intermediate, which polarizes the peroxo bridge and promotes its heterolysis. After O-O bond cleavage, the following reaction steps proceed effortlessly so that the N-nitrosation is accomplished without NO exchange among reaction species.

Interactive Regulation between Aliphatic Hydroxylation and Aromatic Hydroxylation of Thaxtomin D in TxtC: A Theoretical Investigation.

TxtC is an unusual bifunctional cytochrome P450 that is able to perform sequential aliphatic and aromatic hydroxylation of the diketopiperazine substrate thaxtomin D in two distinct sites to produce thaxtomin A. Though the X-ray structure of TxtC complexed with thaxtomin D revealed a binding mode for its aromatic hydroxylation, the preferential hydroxylation site is aliphatic C14. It is thus intriguing to unravel how TxtC accomplishes such two-step catalytic hydroxylation on distinct aliphatic and aromatic carbons and why the aliphatic site is preferred in the hydroxylation step. In this work, by employing molecular docking and molecular dynamics (MD) simulation, we revealed that thaxtomin D could adopt two different conformations in the TxtC active site, which were equal in energy with either the aromatic C20-H or aliphatic C14-H pointing toward the active Cpd I oxyferryl moiety. Further ONIOM calculations indicated that the energy barrier for the rate-limiting hydroxylation step on the aliphatic C14 site was 9.6 kcal/mol more favorable than that on the aromatic C20 site. The hydroxyl group on the monohydroxylated intermediate thaxtomin B C14 site formed hydrogen bonds with Ser280 and Thr385, which induced the l-Phe moiety to rotate around the Cß-Cγ bond of the 4-nitrotryptophan moiety. Thus, it adopted an energetically favorable conformation with aromatic C20 adjacent to the oxyferryl moiety. In addition, the hydroxyl group induced solvent water molecules to enter the active site, which propelled thaxtomin B toward the heme plane and resulted in heme distortion. Based on this geometrical layout, the rate-limiting aromatic hydroxylation energy barrier decreased to 15.4 kcal/mol, which was comparable to that of the thaxtomin D aliphatic hydroxylation process. Our calculations indicated that heme distortion lowered the energy level of the lowest Cpd I α-vacant orbital, which promoted electron transfer in the rate-limiting thaxtomin B aromatic hydroxylation step in TxtC.

Theoretical study on the catalytic mechanism of human deoxyhypusine hydroxylase.

Deoxyhypusine hydroxylase is a critical enzyme for hypusination of eukaryotic translation initiation factor 5A (eIF5A). Human deoxyhypusine hydroxylase (hDOHH) has a nonheme diiron active site that resembles both in structure and function of those found in methane and toluene monooxygenases, bacterial and mammalian ribonucleotide reductases, and stearoyl acyl carrier protein Δ9-desaturase from plants. However, the detailed catalytic mechanism of hDOHH is still unclear. In this work, extensive DFT calculations reveal that the catalytic mechanism of hDOHH consists of four consecutive steps: (1) peroxo isomerization triggered by substrate binding; (2) rate-determining O-O bond cleavage and formation of the [FeIV2(µ-O)2]4+ compound; (3) H atom abstraction from the substrate; and (4) OH rebound to the substrate. This work not only rationalizes the exceptional stability of the diiron(iii)-peroxo complex in hDOHH, but also confirms that hDOHH uses a diamond shape [FeIV2(µ-O)2]4+ core to complete crucial H atom abstraction from the substrate. Our DFT calculations exclude the reaction pathway of hDOHH to use diiron(iii)-peroxo species to directly react with the substrate.

Ruffling drives coproheme decarboxylation by facilitating PCET: a theoretical investigation of ChdC.

Coproheme decarboxylase (ChdC) is an essential enzyme in the coproporphyrin-dependent heme synthesis pathway, which catalyzes oxidative decarboxylation of coproheme at the positions p2 and p4 to generate heme b under the action of hydrogen peroxide. A mysterious characteristic of catalytic mechanism of ChdC is that both of the two decarboxylation sites are located remotely from the iron center of coproheme, which binds with hydrogen peroxide. By using density functional theory calculations, we have studied the coproheme decarboxylation mechanism of ChdC in detail. The calculation results show that in the first step of the catalytic reaction, H2O2 homolysis takes place synergistically with the proton coupled electron transfer process of a tyrosine (Tyr145) residing near p2 propionate. The produced reactive Tyr radical then abstracts a hydrogen atom from the ß carbon of the p2 propionate side chain, which is the rate-limiting step of the whole reaction with a 19.16 kcal mol-1 energy barrier. Finally, through intramolecular electron and proton rearrangement of coproporphyrin, decarboxylation of p2 propionate is accomplished. Our study revealed that the ruffled conformation of coproheme in ChdC is an important structural factor, which facilitates the decarboxylation reaction. We also found that the hydrogen bond chain located below the coproheme ring plays a role to regulate the PCET process of Tyr145. In addition, molecular dynamics simulations discovered that Lys149 is responsible for stabilizing the harderoheme III and positioning the second decarboxylation site p4 to the catalytic Tyr145 site in the decarboxylation reaction of the p4 site.

ONIOM investigations of the heme degradation mechanism by MhuD: the critical function of heme ruffling.

The oxygen-dependent heme utilization degrading enzyme in Mycobacterium tuberculosis (MhuD) uniquely integrates monooxygenase and dioxygenase functions in a single active site. It cannot convert heme to biliverdin as canonical heme oxygenases but generates mycobilin without releasing carbon monoxide. Herein, by employing ONIOM calculations, we investigated the heme degradation mechanism of MhuD. Our calculations revealed that MhuD firstly follows a canonical monooxygenation mechanism to hydroxylate heme on the δ-meso carbon guided by the asparagine residue Asn7, which experiences a 21.2 kcal mol-1 energy barrier in the O-O cleavage rate-limiting step during the conversion process from ferric heme-hydroperoxy species to mycobilin. In the second degradation step, the ruffled conformation of oxoheme (oxoheme is the ferrous π radical complex formed by hydroxyheme experiencing deprotonation in the hydroxyl group and intramolecular electron transfer) imposed by the hydrophobic environment of the enzyme not only inhibits the continuing conversion of oxoheme to biliverdin but also endows the meso-carbons with radical characteristics, which turns the second degradation step to a dioxygenation reaction with 20.4 kcal mol-1 energy barrier. We further analysed the electronic structure change along the reaction process. Our calculation discovered that the ruffled structure of oxoheme is critical to the regiospecificity and even atom location selectivity, as well as the reaction mechanism of the degradation process.

DNA Structural Distortions Induced by a Monofunctional Trinuclear Platinum Complex with Various Cross-Links Using Molecular Dynamics Simulation.

The monofunctional trinuclear platinum complex (MTPC), as a promising antitumor agent, can form MTPC-DNA adducts via bifunctional and trifunctional cross-links. Molecular dynamics simulations were used to investigate DNA structural distortions of the MTPC-DNA adducts. MTPC coordinating to DNA results in the decrease of base-pair thermal stability and DNA structural distortions. It is found that there are more significant DNA structural distortions in the trifunctional cross-link than in the bifunctional cross-link, in the 1,4-GG than in the 1,3-GG cross-link, and in the intrastrand than in the interstrand cross-link with the same spans. The results provide a better understanding of DNA structural distortions induced by MTPC with various cross-links at the nucleotide level and are helpful for exploring novel Pt-based anticancer drugs.

Discovery of Novel Biased Opioid Receptor Ligands through Structure-Based Pharmacophore Virtual Screening and Experiment.

Gi -protein-biased agonists with minimal ß-arrestin recruitment represent opportunities to overcome the serious adverse effects of human mu opioid receptor (µ-OR) agonists and developing alternative and safe treatments for pain. In order to discover novel non-morphinan opioid receptor agonists, we applied hierarchical virtual screening of our in-house database against a pharmacophore based on modeling the active conformations of opioid receptors. We discovered an initial hit compound, a novel µ-OR agonist with a pyrazoloisoquinoline scaffold. We applied computational R-group screening to this compound and synthesized 14 derivatives predicted to be the best. Of these, a new Gi -protein-biased compound, 1-{5-(3-chlorophenyl)-7,8-dimethoxy-3-[4-(methylsulfonyl)benzyl]-3H-pyrazolo[3,4-c]isoquinolin-1-yl}-N,N-dimethylmethanamine, showed an EC50 value of 179 nm against the µ-OR. This resulted in significant pain relief for mice in the phase II period of formalin response tests. This study provides a new strategy to identify diverse sets of promising compounds that might prove useful for the development of drugs that target other G-protein-coupled receptors.

Theoretical study of the mechanism of the manganese catalase KatB.

The mechanism of the H2O2 disproportionation catalyzed by the manganese catalase (MnCat) KatB was studied using the hybrid density functional theory B3LYP and the quantum chemical cluster approach. Compared to the previous mechanistic study at the molecular level for the Thermus thermophilus MnCat (TTC), more modern methodology was used and larger models of increasing sizes were employed with the help of the high-resolution X-ray structure. In the reaction pathway suggested for KatB using the Large chemical model, the O-O homolysis of the first substrate H2O2 occurs through a µ-η1:η1 coordination mode and requires a barrier of 10.9 kcal/mol. In the intermediate state of the bond cleavage, two hydroxides form as terminal ligands of the dimanganese cluster at the Mn2(III,III) oxidation state. One of the two Mn(III)-OH- moieties and a second-sphere tyrosine stabilize the second substrate H2O2 in the second-sphere of the active site via hydrogen bonding interactions. The H2O2, unbound to the metals, is first oxidized into HO2· through a proton-coupled electron transfer (PCET) step with a barrier of 9.5 kcal/mol. After the system switches to the triplet surface, the uncoordinated HO2· replaces the product water terminally bound to the Mn(II) and is then oxidized into O2 spontaneously. Transition states with structural similarities to those obtained for TTC, where µ-η2-OH-/O2- groups play important roles, were found to be higher in energy.

Structural Studies of the 3',3'-cGAMP Riboswitch Induced by Cognate and Noncognate Ligands Using Molecular Dynamics Simulation.

Riboswtich RNAs can control gene expression through the structural change induced by the corresponding small-molecule ligands. Molecular dynamics simulations and free energy calculations on the aptamer domain of the 3',3'-cGAMP riboswitch in the ligand-free, cognate-bound and noncognate-bound states were performed to investigate the structural features of the 3',3'-cGAMP riboswitch induced by the 3',3'-cGAMP ligand and the specificity of ligand recognition. The results revealed that the aptamer of the 3',3'-cGAMP riboswitch in the ligand-free state has a smaller binding pocket and a relatively compact structure versus that in the 3',3'-cGAMP-bound state. The binding of the 3',3'-cGAMP molecule to the 3',3'-cGAMP riboswitch induces the rotation of P1 helix through the allosteric communication from the binding sites pocket containing the J1/2, J1/3 and J2/3 junction to the P1 helix. Simultaneously, these simulations also revealed that the preferential binding of the 3',3'-cGAMP riboswitch to its cognate ligand, 3',3'-cGAMP, over its noncognate ligand, c-di-GMP and c-di-AMP. The J1/2 junction in the 3',3'-cGAMP riboswitch contributing to the specificity of ligand recognition have also been found.

Function Coupling Mechanism of PhuS and HemO in Heme Degradation.

Most bacteria possess only one heme-degrading enzyme for obtaining iron, however few bacteria such as Pseudomonas aeruginosa express two, namely PhuS and HemO. While HemO is a well-known heme oxygenase, previously we discovered that PhuS also possesses heme degradation activity and generates verdoheme, an intermediate of heme breakdown. To understand the coexistence of these two enzymes, using the DFT calculation we reveal that PhuS effectively enhances heme degradation through its participation in heme hydroxylation, the rate limiting reaction. Heme is converted to verdoheme in this reaction and the energy barrier for PhuS is substantially lower than for HemO. Thus, HemO is mainly involved in the ring opening reaction which converts verdoheme to biliverdin and free iron. Our kinetics experiments show that, in the presence of both PhuS and HemO, complete degradation of heme to biliverdin is enhanced. We further show that PhuS is more active than HemO using heme as a substrate and generates more CO. Combined experimental and theoretical results directly identify function coupling of this two-enzyme system, resulting in more efficient heme breakdown and utilization.

Distal Proton Shuttle Mechanism of Ribosome Catalysed Peptide Bond Formation-A Theoretical Study.

In this work, we have investigated a novel distal proton shuttle mechanism of ribosome catalyzed peptide bond formation reaction. The reaction was found to follow a two-step mechanism. A distal water molecule located about 6.1 Šaway from the attacking amine plays as a proton acceptor and results in a charge-separated intermediate that is stabilized by the N terminus of L27 and the A-site A76 5'-phosphate. The ribose A2451 bridges the proton shuttle pathway, thus plays critical role in the reaction. The calculated 27.64 kcal•mol-1 free energy barrier of the distal proton shuttle mechanism is lower than that of eight-membered ring transition state. The distal proton shuttle mechanism studied in this work can provide new insights into the important biological peptide synthesis process.

A theoretical study on oxidative cleavage of olefins to carbonyls catalysed by Fe(iii)-PyBisulidine.

An environmentally friendly new protocol for the selective aerobic cleavage of styrene to carbonyl compounds using the Fe(iii)-PyBisulidine catalyst has been reported recently. The catalyst features several unusual characteristics, such as its high efficiency lies on the ferric center instead of ferrous used by most iron-containing oxygenases and the catalyst specifically oxidizes phenyl-substituted olefins but exhibits no activity on nonconjugated olefins. Herein, we have investigated the mechanism of the oxidative cleavage reaction catalyzed by Fe(iii)-PyBisulidine at the quantum chemistry level. Our computational study shows that the catalyst uses a dioxygen ligation mechanism to activate dioxygen to receive one electron from olefin, which triggers the oxidative cleavage reaction. Our study rationalizes that the Fe(ii)-PyBisulidine catalyst is inactivated because ferrous is unable to raise the oxidizing ability of dioxygen. The exclusive oxidative cleavage of the phenyl-substituted olefin mainly results from the stability of the carbon cation, the orbital symmetry between the conjugated olefin and dioxygen, as well as a lower energy level of HOMO in conjugated olefin.

Influence of Secondary-Structure Folding on the Mutually Exclusive Folding Process of GL5/I27 Protein: Evidence from Molecular Dynamics Simulations.

Mutually exclusive folding proteins are a class of multidomain proteins in which the host domain remains folded while the guest domain is unfolded, and both domains achieve exchange of their folding status by a mutual exclusive folding (MEF) process. We carried out conventional and targeted molecular dynamics simulations for the mutually exclusive folding protein of GL5/I27 to address the MEF transition mechanisms. We constructed two starting models and two targeted models, i.e., the starting models GL5/I27-S and GL5/I27-ST in which the first model involves the host domain GL5 and the secondary-structure unfolded guest domain I27-S, while the second model involves the host domain GL5 and the secondary/tertiary-structure extending guest domain I27-ST, and the target models GL5-S/I27 and GL5-ST/I27 in which GL5-S and GL5-ST represent the secondary-structure unfolding and the secondary/tertiary-structure extending, respectively. We investigated four MEF transition processes from both starting models to both target models. Based on structural changes and the variations of the radius of gyration (Rg) and the fractions of native contacts (Q), the formation of the secondary structure of the I27-guest domain induces significant extending of the GL5-host domain; but the primary shrinking of the tertiary structure of the I27-guest domain causes insignificant extending of the GL5-host domain during the processes. The results indicate that only formation of the secondary structure in the I27-guest domain provides the main driving force for the mutually exclusive folding/unfolding between the I27-guest and GL5-host domains. A special structure as an intermediate with both host and guest domains being folded at the same time was found, which was suggested by the experiment. The analysis of hydrogen bonds and correlation motions supported the studied transition mechanism with the dynamical "tug-of-war" phenomenon.